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ABSTRACT 

HH 666 is an externally irradiated protostellar outflow in the Carina nebula for which 
we present new near-IR [Fe n] spectra obtained with the FIRE spectrograph at Magel¬ 
lan Observatory. Earlier Ha and near-IR [Fe n] imaging revealed that the two emission 
lines trace substantially different morphologies in the inner ~ 40" of the outflow. Ha 
traces a broad cocoon that surrounds the collimated [Fe n] jet that extends through¬ 
out the parent dust pillar. New spectra show that this discrepancy extends to their 
kinematics. Near-IR [Fe n] emission traces steady, fast velocities of ±200 km s -1 from 
the eastern and western limbs of the jet. We compare this to a previously published 
Ha spectrum that reveals a Hubble-flow velocity structure near the jet-driving source. 
New, second-epoch HST/ACS Ha images reveal the lateral spreading of the Ha out¬ 
flow lobe away from the jet axis. Ha proper motions also indicate a sudden increase 
in the mass-loss rate ^ 1000 yr ago, while steady [Fe n] emission throughout the inner 
jet suggest that the burst is ongoing. An accretion burst sustained for ^ 1000 yr is 
an order of magnitude longer than expected for FU Orionis outbursts, but represents 
only a small fraction of the total age of the HH 666 outflow. Altogether, available 
data suggests that [Fe n] traces the highly collimated protostellar jet while Ha traces 
the entrained and irradiated outflow. HH 666 appears to be a missing link between 
bare jets seen in H II regions and entrained molecular outflows seen from embedded 
protostars in more quiescent regions. 
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1 INTRODUCTION 


Over a wide range of initial masses, stars appear to drive 
an outflow at some stage in their formation. Mounting 
evidence calls into question whether there is an upper 
limit on the mass of protostars that may drive colli¬ 
mated jets, despite uncertainty about whether the dom¬ 
inant physics of for mation may c hange b etween low and 
high m asses (e.g. iGarav et a l. 2003; Cara tti o Garatti et al.l 
2008 1: I Guzman et al.l 20121 : iDuarte-Cabral et al.l l20ljj ; 
Reiter V Smithl 120131 ) . These jets will interact with 
and shape - and be shaped by - the circumstellar 
molecular envelope and larger natal cloud. Various jet 
properties and morphologies, such as asym metries (e.g. 
iReipurth et al. 19981: iBa llv V R eipurthll200ll ). wiggles and 
gaps (e.g. iDe Gouveia et al.lll996 ). the degree of collima- 


tion jMundt et al.ll 199 lh , and decreasing knot velocities have 
been explained by interaction with different kinds of envi¬ 
ronments. 


* E-mail: mreiter@as.arizona.edu (MR) 


Interaction between the protostellar jet and ambient gas 
in the star forming cloud has been proposed as a mech¬ 
anism to create molecular outflows. Ra ga V Cabritl (p.993) 
offer a simple analytical model wherein molecular outflows 
are the product of environmental gas being entrained by 
a protostellar jet. While some molec ular outflows are well- 
described by jet-entrainmen t (see, e.g. lZinnecker etahl [l998: 
lArce V Sargentll2005l . I 2 OO 6 I ) , others show a conical morphol¬ 
ogy and broader range of velocities in the outflowing gas that 
have been interpreted as a wide-angle wind launched from 
the protostellar disc (e. g. iNagar et al.lIl997l : iKlaassen et al.l 
120131 : ISalvk et"ai1120141 ) . A full evolutionary scenario may 
involve a combination of these two, with diff erent co mpo¬ 
nents dominating at different times (e.g. IShang et al.l 12006 : 
IZapata et alfe oM). In either case, the outflow will clear ma¬ 
terial along the jet axis, creating a polar cavity as it travels 
through the protostellar envelope and into the surrounding 
molecular cloud. Clearing of the protostellar envelope may 
explain the apparent widening of m olecular o utfl ows a ssoci- 
ated with more evolved sources (|Arce V Sargentl 1200 6). 

In the unobscured environment outside the natal cloud, 
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Table 1 . Observations 


Instrument 

Filter / Position 

Date 

Int. time 

Comment 

Ref. 

ACS/HST 

F658N 

2005 Mar 30 

1000 s 

Ha + [N 11 ] 

Smith et al. (2010i 

ACS /HST 

F658N 

2014 Feb 21 

1000 s 

Ha + [N 11 ] 

this work 

WFC3-UVIS / HS T 

F656N 

2009 Jul 24 

7920 s 

Ha 

Reiter & Smith (2013i 

WFC3-UVIS /HST 

F673N 

2009 Jul 24 

9600 s 

[S 11 ] A6717 + A6731 

Reiter & Smith (2013i 

WFC3-IR/HST 

F126N 

2009 Jul 25 

600 s 

[Fe 11 ] A12567 

Reiter & Smith (2013i 

WFC3-IR/HST 

F164N 

2009 Jul 25 

700 s 

[Fe 11 ] A16435 

Reiter V Smith f20131 

EMMI/JVTT 

Ha 

2003 Mar 11 

2400 s 

P.A. = 293.5° 

Smith et al. (2004i 

FIRE/ Magellan 

HH 666 M 

2013 Feb 22 

600 s 

P.A. = 110° 

this work 

FIRE /Magellan 

HH 666 IRS 

2013 Feb 22 

600 s 

P.A. = 110° 

this work 

FIRE/ Magellan 

HH 666 O 

2013 Feb 22 

600 s 

P.A. = 115° 

this work 


optical and near-IR emissio n lines can be u sed as a tracer 
of the protostellar jet (e.g. ISeoulveda et al.ll201lh . Longer 
wavelength emission lines from shock-excited molecules, 
on the other hand, sample embedded parts of the en- 
tr ained outflow propag atin g within a molecular cloud (e.g. 
iNoriega-Cresoo et al.ll2004l ). In this paper, we refer to the 
jet as the fast, highly collimated stream of gas launched 
near the poles of the protostar and the outflow as the 
slower, wide-angle component of outflowing gas that ei¬ 
ther originated in the protostellar disc, or in the ambient 
medium that was entrained by shocks in the passing jet. 


(e.g. Beltran et al.1 

2008 

Beuther et al.| 20081; Forbrich et al. 

2009: Fuente et al.1 

2009 

Takahaski V Ho 20l3i. The under- 


lying collimated jets are primarily atomic and remain in¬ 
visible unless excited by shocks or external irradiation. For 
both excitation mechanisms, bare jets emerging from more 
evolved, unembed ded sources will emit at visual and UV 
wave l engths (e.g. Reipurth et Ml Jl99S; iBallv &; Reipurthl 


I 2 OOII : iMcGroartv et al.ll2004| j. Observing both the jet and 
the outflow in a single source is difficult, but understanding 
how the jet interacts with the surrounding cloud is impor¬ 
tant for constraining local momentum feedback. Partially 
embedded jets offer one way to bridge these two regimes. For 
example, a protostar residing near the edge of a cloud or in 
a small, isolated globule may drive a jet into the relatively 
diffuse environment outside the cloud. At the same time, the 
counter-jet plows deeper into the molecular cloud where its 
presence must be inferred from the shape of the outflow lobe 
it creates. Jets and outflow lobes have been seen to coex- 


ist in several protostellar sources fe.g. HH lll:lNagar et al 


2004: lArce et al, 


Arce &; Goodman l2001al j. 


1997: iLeFloch et ahll2007l : HH 46/47: [Noriega- Cresp o et al 


2013; and HH 300: IReipurth et al.l l2000l : 


In the Carina nebula, HH 666 emerges from a translu¬ 
cent pillar into a giant H.II region created by more th an ~ 65 
O-type stars (see Smithl l2006ah . ISmith et al.l (|2004| j discov¬ 
ered the outflow in ground-based narrowband optical images 
and identified several large knots that extend more than 4 
pc on the sky (features HH 666 D, A, E, M, O, N, I, and 
C). T he mass-loss rates (ISmith e t al. 2010 : Reiter fe Smithl 
I2013IJ and outflow velocities (|Reiter fe Smithl l2014h mea- 
sured in HH 666 are among the highest of the jets in Carina, 
making it one of the most power ful jets known (see, e.g. Fig¬ 
ure 14 in lEllerbroek et aLll2013l h Harsh UV radiation from 
the many massive stars in Carina permeates the region, il¬ 


luminating the dust pillar from which HH 666 emerges and 
irradiating the body of the jet itself. External irradiation 
lights up neutral material in the jet that is not excited in 
shocks, revealing portions of the jet that would otherwise re¬ 
main invisible. This allows the physical properties of the jet 
to be derived using the diagnostics of photoionized gas (e.g. 
IBallv et al.ll2006l : ISmith et al.ll2010h . However, massive star 
feedback illuminates indiscriminately, irradiating both the 
jet and the surrounding environment, which can complicate 
the interpretation of the observed emission. Non-uniform 
density structures in a photoionized jet can further compli¬ 
cate the interpretation. Fortunately, multi-wavelength ob¬ 
servations allow us to disentangle the jet from its complex 
environment. 


ISmith et al.1 (|2004l i and lReiter V Smithl (|2013h showed 
that near-IR [Fe 11 ] emission from HH 666 traces the jet in¬ 
side the ~ 40" wide pillar, connecting the large-scale jet seen 
in Ha to the IR-bright protostar that drives it (HH 666 IRS). 
High angular resolution near-IR narrowband [Fe 11 ] images 
obtained with Wide-Field Camera 3 (WFC3-IR) on the Hub¬ 
ble Space Telescope (HST) rev eal additional morpho logical 
complexities (see Figure fTI and lReiter fe Smithl I2013T ). Near 
the driving source, [Fe 11 ] emission traces portions of the jet 
that are not seen in Ha. In the western limb of the inner jet, 
designated HH 666 M, [Fe 11 ] emission bisects a dark cavity 
to the northwest of HH 666 IRS, clearly connecting the jet 
to its driving source. [Fe 11 ] emission from HH 666 O, the 
eastern limb of the inner jet, traces a collimated jet body 
that appears to be surrounded by_a broad cocoon of Ha 
emission. Curiously, [Smith et al.l (2004) see a Hubble-like 
velocity structure in Ha and [S 11 ] spectra of the inner parts 
of HH 666 (features HH 666 M and O). Only two other pro¬ 
tostellar jets have shown Hu bble- li ke veloci ty structures in 
optical spectra (H H 83. IReipurth et all Il989l : and HH 900, 
iReiter et alJl2013 ). 


In this paper, we present new ground-based near-IR 
[Fe 11 ] spectroscopy and a second epoch of HST /ACS Ha 
images. When combined with extant data, the HST images 
allow us to compare the kinematics of these two morpholog¬ 
ically distinct components. Together, these new data offer 
a self-consistent explanation for the curious morphological 
and kinematic differences between the Ha and [Fe 11 ] emis¬ 
sion from HH 666. 
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Figure 1. Top: New HST /ACS Ha image of HH 666 with labels and arrows pointing o ut features discussed in the text and used to 
measure proper motions (a). Middle: H5T/WFC3-UVIS H a (b) and [S il] (c) images from lReiter fc Smithl (l2013li with an unsharp mask 
applied to improve the contrast between faint jet features and bright emission from the surface of the dust pillar. Narrowband, near-IR 
[Fe il] images with an unsharp mask applied to improve contrast are shown in the row below (panels d and e). Bottom: An H<a-[S il] 
difference image (f) showing the offset between the two emission lines (excess Ha emission appears white, [S il] emission is black). Small 
images on either side show a zoomed-in view of the two places in the innermost jet where the two emission lines are significantly offset. 
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2 OBSERVATIONS 

2.1 Near-IR spectroscopy 

Near-IR spectra of HH 666 were obtained with the Folded- 
Port Infrared Echellette spectrograph (FIRE, ISimcoe et al.l 
teoiah on the Baade 6.5-m Magellan telescope on 22 February 
2013. Conditions were photometric with thin clouds and the 
seeing was ~ 0'/6. FIRE captures medium resolution spec¬ 
tra from 0.8 — 2.5 pm simultaneously. The echelle slit is 7" 
long and our data were obtained with a slit width of 0'/6. 
This corresponds to a velocity resolution of 50 km s -1 . Un¬ 
certainties in the measured radial velocities are shown in 
Figure [3] We used five slit positions to target the brightest 
[Fe n] emission in the inner ~ 40^ of the HH 666 jet and its 
driving source. Slit positions are shown in Figure [3] 

We obtained separate frames offset ~ 5" from the jet 
axis for the four pointings along HH 666 M and HH 666 O to 
subtract sky lines and emission from the H n region. For the 
central slit position that contains the driving source (see Fig¬ 
ure n, we employ an ABBA nod putting the driving source 
at opposite ends of the slit and remove sky emission by sub¬ 
tracting complementary nods. The total integration time for 
each slit position is 600 s. The data were reduced using the 
firehose IDL pipeline. This data reduction package pro¬ 
vides flat-fielding, sky-subtraction, and wavelength calibra¬ 
tion of the 2D FIRE spectra. Wavelength calibration was 
done using a ThAr lamp. We extract the [Fe II] position- 
velocity diagrams from the reduced, calibrated 2D spectra. 
The reported velocities are in the heliocentric velocity frame. 


2.2 New Ha images 

We present new Ha images obtained with HST/ACS on 
21 February 2 014. These images d uplicate the observational 
setup used bv lSmith et al . (12010 0 to observe HH 666 on 30 
March 2005, providing two epochs with a time baseline of 
8.9 yr (see Table 1). Both epochs were observed using the 
F658N filter (which transmits both Ha and [N II] A6583) 
with a total exposure time of 1000 s. A three-p osition linear 
offset pattern was used to fill ACS chip gaps. See lSmith et al .1 
( 2010 0 for a more complete description of the observational 
setup. 

By replicating the observational setup of the first- 
epoch images, we reduce systematic uncertainties in the 
image registration and measurement of individual blob 
proper motions. Our image regi stration procedure is si mi¬ 
lar to the method described in lAnderso iL_eUa]J ( 2008a, b), 
lAnderson &; van der MareJ (|2010h , and ISohn et al.l (|20ld ) . 
This method uses PSF photometry of pipeline-calibrated im¬ 
ages that have been bias-subtracted, flat-fielded, and CTE- 
corrected, but have not been resampled to correct for ge¬ 
ometric distortion (denoted with the file extension flc). 
These images are better suited for high-accuracy PSF fitting 
than drizzled images (denoted by the file extension drc), and 
are therefore preferable to measure precise relative proper 
motions from HST data. 

We construct an initial reference frame from the drc 
images using the astrometric information in their headers. 
Basic centroid positions of bright, relatively-isolated stars in 
the drc images are compiled into a reference frame that has a 
50 mas pixel -1 scale and is aligned with the y axis pointing 


north, like that of lAnderson et all ( 2008a f. We then per¬ 
form PSF photometry on e ach flc exposure using the pro¬ 
gram img2xym_WFC. 09x10 (Anderson &; Kingl 120061 ). which 
employs a spatially-variable library of PSFs. Stellar posi¬ 
tions are corrected afterward for geometri c dist ortion using 
the distortion correction from [Xndersonl ( 2006 ). Next, six- 
parameter linear transformations are computed from each 
distortion-corrected flc frame to the initial reference frame 
via the positions of the stars found in both. For improved ac¬ 
curacy, we recalculate the transformation after replacing the 
initial reference-frame positions of the stars with their aver¬ 
age PSF-based positions transformed from the flc frames. 

Finally, we construct a stacked image for each epoch by 
resampling the flc images into the reference frame, allowing 
us to measure the motions of extended sou rces. The details 
of the stacking algorithm can be found in lAnderson et all 
(2008a). As both epochs are transformed to the same ref¬ 
erence frame, object positions are directly comparable. The 
image alignment precision is primarily limited by the inter¬ 
nal accuracy of the ref erence frame and distortion correc¬ 
tion 0.015 pixels, see lAnderson fe Kindl2006l : ISohn et al.l 
120121 ). which over our time baseline corresponds to 0.084 
mas yr - 1 or 0.92 km s -1 at the distance of Carina (2.3 kpc, 
ISmithll2006bl ). 


To measure proper motions of bright knots in the in¬ 
ner jet, we use our implemen tation o_f_the modified cross¬ 
correlation techn i ciue used by Currie et al.l (Il996|. see also 
Morse et al 200lUHartigan et al 200ll . and lReiter fc Smithl 


2014) . First. we correct for local diffuse H 11 region emis¬ 


sion by subtracting a median-filtered image from each epoch. 
Next, we extract bright jet knots using a box size that is op¬ 
timized for each feature (see Figure H|). To determine the 
pixel offset between the two epochs, we generate an array 
that contains the total of the square of the difference be¬ 
tween two images for each shift of a small box containing 
the jet feature relative to a reference image. The best match 
between the two images will produce the smallest value in 
the array. 

With the identical observational setup and longer time 
baseline between the two epochs of HST/ACS data, we 
achieve finer vel ocity resolution than that obtained in our 
previous study ([Reiter &; Smithl 120141 ) using HST/ WFC3- 
UVIS images taken 4.32 yr after the first-epoch HST /ACS 
images. Using the same camera, filter, and approximate lo¬ 
cation of the source on the detector (and thus same dis¬ 
tortion correction) minimizes systematic uncertainties and 
allows us to measure proper motions of bright nebular jet 
features with velocities greater than ~ 10 km s -1 . To de¬ 
termine the uncertainty in the pixel offset between the two 
epochs, we explore the proper motion obtained using differ¬ 
ent box sizes around each jet feature. We also determine the 
maximum offset between the two images that still results 
in good alignment, identified by residuals in a subtraction 
image on the order of the noise. From this, we determine 
an average proper motion uncertainty of ~ 10 km s -1 , al¬ 
though this depends on the signal-to-noise of the feature (see 
Table 2). 


2.3 Previously published images and spectroscopy 

For comparison to the new Ha image and [Fe II] spectra, we 
include HST/WFC3 narrowband Ha, [S 11], and [Fe 11] 1.26 
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p m and 1.64 pm images . These were discussed previously 
bv lReiter fc Smith ( 2013 ). Full observational details can be 
found in that paper. 

W e also inclu de the Ha spectrum of HH 666 from 
ISmith et al.l ([20041 ). Longslit Ha observations were obtained 
over a slit length of ~ 6b With the cross-dispersing grism 
and a 1" slit width, the H a spectrum has a velocity resolu¬ 
tion of ~ 10 km s _1 . More com plet e d etails of the observa¬ 
tion and reduction are given bv ISmith et alJ (120041 ). 


3 RESULTS 

3.1 Morphology in images 

High angular resolution images from HST reveal the com¬ 
plex morphology of HH 666, particularly in the inner ~ 40 /7 
of the jet. Ho emission from HH 666 M begins ~ 10 7/ to the 
northwest of the driving source tracing a broad arch around 
a possible evacuated cavity (see Figure [T]) before quickly fo¬ 
cusing into a narrow beam. On the opposite side of the driv¬ 
ing source (to the southeast), a few Ha-bright knots trace 
the inner jet. Two bright parabolic arcs with extended wings 
reach back toward the driving source creating a broad co¬ 
coon morphology (features 02 and 03 in Figure S}. Apply¬ 
ing an unsharp mask to the Ha image reveals collimated jet 
emission just inside the terminus of HH 666 O (inside feature 
03, see Figures [T| and H]). Comparing Ha and [S n] images 
from HST /WFC3-UVIS demonstrates that the morphology 
of these two l ines is remarkably sim ilar in externally irradi¬ 
ated HH jets ([Reiter fe Smithl[2013l ). In purely shock-excited 
gas, however, the locat ion and relat i ve intensity of the two 
lines will be offset (lHartigan et al.l Il994l : iHeathcote et al.l 
1996). Thus, in the harsh UV environment of the Carina 
nebula, photoionization dominates the excitation of the jet, 
leading to significantly brighter Ha emission compared to 
[S II], with the typical ratio in HH 666 of [S n]/Ha ~ 0.03. 

There are only two places in the jet where bright [S n] 
emission diverges from the morphology of the Ha emission 
(see Figure [T] bottom panel). [S n] emission extends inside 
the first arc of Ha emission from the jet on either side of 
HH 666 IRS. To the northwest (HH 666 M), [S n] emission 
falls just inside the arc of Ha emission that appears to trace 
a cavity, extending ~ 0 7/ 25 behind the Ha. On the other side 
of the driving source in HH 666 O, [S il] is much more pro¬ 
nounced in a difference image, extending nearly 1" behind 
the Ha bright knot that marks the emergence of HH 666 O. 
In a post-shock cooling zone, bright [S n] emission will 
trail an Ha-bright arc, with the specific geometry depending 
on the density and ionization fractio n of the preshock gas , 
and the velocity of t he shock itself (jHartigan et alJ 1 1994 
IHeathcote et alJ fl996) . Such offsets have been observed in 
high-spatial resolution images of nearby HH jets, and inter¬ 
preted as spatially resolved post-shock cooling zones (e.g. 
Heat hcote et al.|[l996l : iReipurth et alJfl997al : lHartigan et al.l 

20111 ). Therefore, interpreting HH 666 as a purely photoion- 
ized jet would be an oversimplification. 

Near-IR [Fe n] emission from HH 666 tells a remark¬ 
ably differ ent story from that told by the optical emis¬ 
sion ([Reiter &; Smithl l2013h . A stream of [Fe n] emission 
in HH 666 M bisects the dark cavity seen in Ha images, 
connecting the jet to the driving source. Outside the Ha- 
bright arc bordering the cavity (features M8, M5a, M10, 




offset from YSO [arcsec] 

Figure 2. (a) [Fe il] 1.26 pm image of HH 666 with boxes showing 
the tracing through the jet to measure the [Fe il] ratio 7 Z plotted 
below, (b) Plot of the [Fe il] 1.26 pm / 1.64 pm ratio 7Z along the 
inner jet. The shaded area surrounding the ratio line indicates the 
standard deviation, (c) Tracing of the [Fe il] 1.26 pm and [Fe il] 
1.64 pm flux through the same portions of the inner jet. 

M12), [Fe n] emission follows the same morphology as Ha. 
Morphological discrepancies are even more striking on the 
opposite side of the jet in HH 666 O. [Fe n] emission does 
not connect the Ha jet all the way to the driving source, 
instead beginning ~ 6 7/ away from the driving source, closer 
than the ~ 11" offset in the Ha emission. From this point, 
[Fe n] emission, although clumpy, extends continuously to 
the bright shock just outside the pillar edge. There is no Ha 
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or [S n] emission from the collimated jet inside the first Ho 
arc in HH 666 O (feature 02 in Figured)). Outside this point, 
Ho and [S n] emission trace both the cocoon and the jet. Op¬ 
tical emission from the outer jet traces the same morphology 
as the [Fe n] jet that bisects the larger ‘sheath.’ Ho, [S n], 
and [Fe n] emission are all bright at the tip of HH 666 O 
just outside the edge of the parent pillar (features 09, 08, 
03, 04, 05 in Figured]). 


3.2 [Fe II] ratio 

We also use the [Fe n] images from iReiter &; Smithl ( 2013 1 
to measure the flux ratio of the two observed lines, 1Z = 
A12567/A16435. Both lines originate from the same upper 
level, so their intrinsic ratio is determined by atomic physics. 
The observed ratio will decrease as the column of absorbing 
material increases, s o 1Z p rovides a measure of the reddening. 
ISmith &, Hartiganl (|2006h measured 7 Z — 1. 49 in P Cygni 
while more recent measurements of HH 1 by iGiannini et al.l 
(|2015h yielded 1Z — 1.11 — 1.20. Since we do not have simul¬ 
taneous off-line continuum images, we estimate local con¬ 
tinuum emission from adjacent regions in the [Fe n] images 
themselves. We plot 1Z along the length of the inner jet in 
Figure [2) 

The median flux ratio through the jet is 1Z ~ 0.98±0.15. 
This corresponds to an A y of at least ~ 0.85 mag assuming 
7Z — 1.49 and ~ 0.25 — 0.41 mag assuming 7Z — 1.11 
1.20. respectively. Both sides of the jet show an increase 
in the ratio along the jet moving toward the edge of the 
pillar. The flux ratio is lowest near the driving source, with 
a median 1Z ~ 0.63 , corre sponding to Av > 1.75 mag for the 
ISmith &; Hartigan (2006) val ue of 7 Z, or Ay > 1.15 — 1.31 
mag using the IGiannini et al.l (|2015l j values. In contrast, the 
ratio is closest to the intrinsic value where we do not observe 
strong [Fe n] emission from the jet (e.g. outside feature 04 
and between features M13 and M14 in Figure [2}. 


3.3 Velocities from spectra 

In the HH 666 discovery paper, ISmith et al.l (|2004l ) pub¬ 
lished longslit Ho and [S n] spectra that confirmed that 
bright knots to the east and west of the dust pillar are part 
of a coherent bipolar outflow. Doppler velocities in Ho and 
[S n] spectra show that the western limb of the jet (including 
HH 666 M) is blueshifted and the eastern limb (HH 666 O) is 
redshiftech with velocities from the inner jet reaching ±250 
km s _ 1 ([Smith et al.l l2004h . The rapid increase in veloc¬ 
ity away from the driving source on both sides of the jet 
shows an almost Hubble-like velocity structure. In a Hub¬ 
ble flow, ve l ocity is proportional to the distance, v oc d. 
[Sm ith et al. (|2004| j derive a “Hubble constant” for HH 666 
of Hhh 666 — 1200(tanz) km s -1 pc -1 , where i is the incli¬ 
nation for the fast inner jet. A few of the more distant knots 
in HH 666 also show an increase in velocity with distance 
from the driving source. Both Ho and [S n] emission traces 
this velocity structure, mirroring the close match of their 
morphology in images. 

The kinematics traced by [Fe n] emission are as dis¬ 
crepant as the morphology in images. Unlike the Her kine¬ 
matics of HH 666, the [Fe n] emission traces steady velocities 
on either side of the jet. An absence of sharp velocity drops 


that suggests that no strong shocks decelerate the jet. Veloc¬ 
ities traced by near-IR [Fe n] are ±200 km s -1 throughout 
the inner jet (see Figure [3}. 


3.4 New proper motions 

We compare new ACS Ho images obtained in 2 014 to the 
firstj epoch images obtained ~ 9 yrs earlier bv ISmith et al .1 
( 2010 ) and measure proper motions of jet features. In gen¬ 
eral, the measu r ed ve locities agree with those found by 
IReiter &; Smithl (120141 ) within the uncertainties (see Ta¬ 
ble 2), although the better match of the instrument and 
filter used in this work permit higher precision than in our 
previous work. 

With the longer time baseline and better match of the 
observational setup, we can also measure proper motions 
of slower features and off-axis motions with greater fidelity. 
The most remarkable motions we uncover with the latest 
epoch shows clear lateral expansion of the cocoon seen in Ho 
emission (see Figure [4]). In particular, we measure a trans¬ 
verse velocity for feature M12 of 64 km s -1 and it appears to 
move at an angle ~ 10° away from the jet axis defined by the 
motion of features M8, M5a, M10. Similar lateral expansion 
in features 06, 07 from HH 666 O traces the slow expan¬ 
sion of the H a sheath surrounding the [Fe n] jet. A slower 
velocity and an off-axis propagation directio n are predicted 
for the expanding wings of a bow shock (e.g. lOstriker et al.l 

1 200 il l. 


4 RESOLVING MORPHOLOGICAL AND 
KINEMATIC DISCREPANCIES 

Given the striking disparity in the morphology and kinemat¬ 
ics traced by Ha and [Fe n] emission in both images and 
spectra, we propose that these two different emission lines 
trace physically distinct components of the HH 666 outflow. 
We illustrate the anatomy of a dense, externally irradiated 
protostellar outflow in Figure [5) and discuss the evidence for 
this picture below. 


4.1 Redshifted lobe: HH 666 O 

In HH 666 O, the redshifted part of the outflow to the south¬ 
east of the protostar, Ho emission appears to form a cocoon 
that surrounds the [Fe n] jet core. The large spatial offset 
of the Ho emission from the [Fe n] (~ 2") suggests that 
it traces an entirely different component of the outflow. In 
fact, the thin arcuate morphology of the Ho and [S n] emis¬ 
sion along the face of the pillar resembles the molecular out- 


flow lobes seen at longer (mm) waveleng 

;ths ( 

e.g. HH 211, 

Gueth & Guilloteau 1999; Dionatos et al. 

2010 

;[Tappe et al. 

2013). Prompt entrainment of ambient 

material via bow 


shocks from episod ic jets has been invoke d to e xplain nested 
outflow shells (e.g. iNaravanan &; Walker! 1 1996l b Jets in Ca¬ 
rina that are driven by protostars embedded in a dust pillar 
(more than half the jets discovered by ISmith et al.l l2010h 
will entrain material as they exit the cloud. 

Kinematics in the Ho spectrum also support the inter¬ 
pretation that Ho emission traces an externally irradiated 
outflow shell in HH 666. Fast, ionized material at the head 
of the bow shock that creates this outflow shell corresponds 
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Figure 3. (a) H a position-velocity diagram of HH 666 from lSmith et al . (12004 1. (b) New, near-IR [Fe il] 1.26 pm spectra obtained with 
five separate slit positions shown to the right. Alongside the [Fe il] position-velocity diagram, we p lot the median line velocity over steps 
of ~ Of !6 (c). Blue boxes show the Ha velocities for HH 666 M and O determined bv ISmith et al. ( 2004 1. (d) FIRE slit positions shown 
as red boxes on a [Fe il] 1.26pm + 1.64 pm image from HS T /WFC3-IR. The dotted line indicates the slit position used to obtain the 
longslit Ha spectrum. 


to the highest velocities seen in the Ha spectrum. Velocities 
observed closer to the driving source are lower because they 
are dominated by emission from the slower wings of the bow 
shock. The slowest velocities will be observed farthest from 
the apex of the bow shock, and therefore closest to the d riv- 
ing source (see, e.g. Figure 1 in lArce fc Goodman!l2002l l. In 
this picture, the bow shock dominates entrainment of the 
ambient gas and will create the observed Hubble-like veloc¬ 
ity structure, as seen in the po sition-velocity diagrams o f 
molecular outflows (e.g. HH 300. lArce V Goodmanll2001bl l . 
In contrast to the Ha emission, Doppler velocities in the 
collimated [Fe 11 ] jet are remarkably constant throughout 
HH 666 O, remaining close to +200 km s -1 , as expected for 
an unimpeded jet propagating away from its driving source. 
We do not detect evidence for a Hubble-flow in the [Fe 11 ] 
emission from the jet. If environmental interaction is dom¬ 
inated by the first impulse from the jet, portions of the jet 
traced by [Fe 11 ] that travel in the wake of the initial bow 
shock will transfer little momentum to the environment. 

Some optical emission from the jet body in the outer 
portions of HH 666 O, near the leading bow shock (feature 
03 in Figure S] also see Figure [I]) further supports the idea 
that there are two externally irradiated outflow components 
visible in HH 666. Ha, [S 11 ], and [Fe 11 ] emission begin to 
coincide just beyond the edge of the dust pillar, marking 
the first place where the newly liberated jet is laid bare 
to ionizing radiation. Optical emission from the jet beyond 
the pillar edge suggests a decrease in the intervening opti¬ 


cal depth, and the increase of the [Fe 11 ] ratio 1Z toward the 
terminus of HH 666 O supports this interpretation. How¬ 
ever, the ratio only approaches the intrinsic value beyond 
the end of the jet, suggesting that the observed reddening is 
primarily from the jet itself (see Figure [2}. 


4.2 Blueshifted lobe: HH 666 M 

The blueshifted part of the outflow, HH 666 M, that extends 
to the northwest of the protostar reveals a similar picture 
to HH 666 O. Close to the driving source (within < 10"), 
the behavior of the Ha and [Fe 11 ] emission in HH 666 M 
resembles that of HH 666 O. A narrow column of knotty 
[Fe 11 ] emission bisects broad, arc-like Ha emission before 
the morphologies traced by the two lines converge at larger 
distances from the protostar. New Ha images show that the 
head of the Ha arc moves along the jet axis while its wings 
expand laterally away from it. Features M8, M5a, M10 (see 
Figure 2} move primarily along the jet axis with a transverse 
velocity of ~ 150 km s -1 . Feature M12 travels at an ~ 10° 
angle away from the jet axis, indicating that this emission 
feature traces the expansion of the Ha arc away from the jet 
axis. Expanding Ha emission suggests that the arc may be 
the bow shock where HH 666 M breaks out of the front side 
of the parent cloud (approximately at the location of feature 
M9). In this picture, the ‘dark patch’ to the northwest of the 
driving source is a cavity cleared by the jet as it plows out of 
the molecular cloud, leaving a deficit of material behind the 
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Table 2. Proper motions in the HH 666 inner jet 


Object 

(5x 

mas 

<5y 

mas 

VT a 
[km s -1 ] 

ViJ b 
[km s -1 ] 

velocity 0 
[km s" 1 ] 

a 

[degrees] 

age d 

yr 

HH 666 Ml 

-125 (1) 

-34 (1) 

159 (4) 

-104 (6) 

190 (7) 

-33 (2) 

1108 (34) 

HH 666 M14 

-122 (3) 

-30 (2) 

154 (5) 


161 (5) 


1097 (42) 

HH 666 M13 

-122 (1) 

-33 (3) 

155 (4) 


162 (4) 


850 (27) 

HH 666 M2 

-124 (1) 

-40 (1) 

159 (4) 


167 (4) 


789 (25) 

HH 666 M3 

-113 (0.5) 

-32 (0.2) 

143 (3) 

-131 (7) 

194 (8) 

-42 (1) 

779 (24) 

HH 666 M4 

-131 (1) 

-25 (2) 

163 (4) 


170 (4) 


593 (19) 

HH 666 M9 

-126 (2) 

-20 (4) 

156 (4) 


164 (4) 


568 (19) 

HH 666 M8 

-123 (1) 

-41 (1) 

159 (4) 


166 (4) 


537 (17) 

HH 666 M5a 

-117 (1) 

-31 (2) 

149 (3) 

-193 (6) 

243 (7) 

-52 (1) 

576 (18) 

HH 666 M10 

-114 (2) 

-37 (1) 

147 (4) 

-193 (6) 

242 (7) 

-53 (1) 

472 (16) 

HH 666 M7 

-124 (1) 

-37 (1) 

158 (4) 


166 (4) 


405 (13) 

HH 666 Mil 

-95 (3) 

-19 (3) 

118 (5) 


124 (5) 


462 (21) 

HH 666 M12 

-52 (0.1) 

-7 (0.3) 

64 (1) 


67 (2) 


1193 (37) 

HH 666 Ol 

107 (2) 

66 (1) 

154 (4) 


162 (4) 


741 (24) 

HH 666 O10 

123 (3) 

47 (3) 

162 (5) 


169 (6) 


895 (35) 

HH 666 02 

88 (3) 

63 (6) 

132 (6) 

223 (18) 

259 (19) 

59 (4) 

1249 (66) 

HH 666 09 

117 (3) 

94 (5) 

184 (7) 

57 (8) 

192 (10) 

17 (1) 

989 (41) 

HH 666 08 

134 (2) 

67 (2) 

184 (5) 

57 (8) 

192 (9) 

17 (0.4) 

1017 (34) 

HH 666 03 

108 (1) 

54 (1) 

148 (4) 


155 (4) 


1363 (44) 

HH 666 04 

91 (1) 

59 (2) 

133 (3) 

121 (4) 

180 (5) 

42 (1) 

1564 (52) 

HH 666 05 

108 (2) 

47 (2) 

144 (4) 


151 (4) 


1434 (52) 

HH 666 06 

50 (67) 

52 (52) 

88 (73) 


92 (77) 


2156 (1792) 

HH 666 07 

9 (3) 

35 (2) 

45 (3) 


47 (3) 


3210 (194) 


Uncertainties for each quantity are listed in parenthesis. 

a The transverse velocity, assuming a distance of 2.3 kpc, measured over At = 8.9 yr. 
b The radial velocity. 

c The total space velocity, assuming the average inclination of 42° where the radial velocity could not be matched. 
d Time for the object to reach its current position at the measured velocity, assuming ballistic motion. 


bow sh ock. ISm ith et ah (l2Qlo|) a l so suggested that a Canto 
nozzle ([Canto fe Rodriguezl 1980l : ICantd et al.ll 198 lh may be 
responsible for refocusing the jet at this point. 


Spectra support the interpretation that the Ha arc 
traces the tip of the expanding outflow cavity. We observe 
a Hubble wedge in HH 666 M, with Ha emission tracing 
increasing velocities up until a sharp drop corresponding to 
features M5a, M10, M12 (see Figure [3}. Inside feature M9, 
closer to the driving source, there is a clear deficit of emission 
from the globule. However, at the same place along the slit 
where globule emission decreases, we see blueshifted emis¬ 
sion that extends smoothly toward feature M9. Thus, it ap¬ 
pears that H a emission from the inner portion of HH 666 M 
also traces the expansion of the outflow shell. 


As in HH 666 O, Doppler velocities measured in the 
[Fe n] spectra of HH 666 M remain steady, with velocities 
close to —200 km s -1 throughout. Fast, oppositely directed 
[Fe n] emission is what we expect from a spatially resolved 
steady jet. This structure has been seen in othe r jets (e.g. 


Garcia-Lopez et all 120081 : iMelnikov et al.l 120091 : iLiu et all 
2012), inc luding those driven bv inter mediate-mass proto¬ 
stars (e.g. lEllerbroek et al.ll2013 . l2014h . Hubble-like veloci¬ 
ties have been observed in protostellar outflows , but usually 
in IR/mm tracers of molec ular gas that has been ac celerated 
by an underlying jet (e.g. lArce fe Goodmanl Il2001b h Only 
two other jets show Hubble-lik e velocities in Ha emission 
- HH 83 (|Reipurth et al.l [l989h and HH 900 ([Reiter et al.l 

HqH). 


A Hubble-like velocity pattern in HH 666 M further 
suggests that the dark patch to the northwest of the driv¬ 
ing source is an evacuated cavity in the dust pillar cleared 
by the jet. If this is correct, HH 666 O should create a cor¬ 
responding cavity to the southeast of the driving source as 
it propagates deeper into the cloud. To test this, we took 
an intensity tracing through the pillar along the axis of the 
jet to search for a decrease in the ambient Ha emission on 
either side of the driving source (see Figure[6|. Indeed, inten¬ 
sity tracings reveal dips in the Ha intensity along the pillar 
extending ~ 5 /x on either side of the driving source. The 
emission trough is shallower to the southeast of the driving 
source, as expected for more foreground emission from the 
pillar itself. 


Beyond the Ha arc (west of feature M10, outside the 
globule), Ha and [Fe n] emission converges, thereafter trac¬ 
ing the same morphology. The [Fe n] ratio remains relatively 
constant (see Figure [2]), and velocities traced by both lines 
remain relatively steady (no strong shocks). Doppler veloci¬ 
ties in [Fe n] spectra, however, are nearly 100 km s -1 faster 
than those traced by Ha. Low- and high-velocity compo¬ 
nents have been observed in jets from low-mass T Tauri stars 
with velocity differences between the core of the jet and the 
slightly off-axis win gs of as much as ^ 100 km s" 1 in som e 
cases (e.g. DG Tau. iBacciotti et ahll2000l : IPvo et ahll2003h . 
If the shell of material that is traced by Ha inside the glob¬ 
ule is pushed back against the jet by feedback from massive 
stars in the region, then the two components may overlap 
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Figure 4. Top: Difference image of the first and second epochs of HST /ACS imaging (a). Middle: New HST /ACS Ha image of 
HH 666 O (b) and HH 666 M (c). Boxes indicate features used to measure proper motions with the name of the feature and the measured 
transverse velocity labelled alongside. Bottom: Same boxes as the middle panel with vectors to indicate the magnitude and direction of 
the measured motion (d and e). 


spatially in images while the entrained material traced by 
Ha has decelerated as compared to the velocity of the [Fe n] 
jet. 


4.3 An externally irradiated jet and outflow 

For embedded jets that entrain an outflow as they plow 
through a dusty pillar, UV radiation from nearby massive 
stars will illuminate the shell of the outflow in addition to 
the jet. We argue that Ho proper motions and Doppler ve¬ 
locities trace the expansion and propagation of the exter¬ 
nally irradiated outflow lobe entrained by the jet, but not 
the jet core. In this picture, we expect larger proper motions 
from the jet itself (~ 18 mas yr -1 for a velocity of 200 km 
s _ 1 ), although additional [Fe n] imaging at a later epoch is 
required to test this hypothesis. 

Our data show that Ha emission traces entrained out¬ 
flowing gas with the same Hubble-like velocity structure seen 
in molecular outflows associated with o ther HH jets (e.g. 
iLee et al1l200ll : lArce &; Goodmanll2001bh . However, the ex¬ 
treme UV radiation in the Carina nebula illuminates en¬ 


trained gas in HH 666, leading to Ha emission from the out¬ 
flow cocoon. No H 2 emission that is unambigu ously associ¬ 
ated with the jet has been observed in HH 666 ([Smith et al.l 
2004: lHartigan et al.ll2015h . although weak H 2 emission will 
be difficult to discern against the bright background of the 
photodissociation region (PDR). 


5 AN EXPANDING OUTFLOW SHELL FROM 

AN ONGOING JET BURST 

Both morphology and kinematics reveal an Ha outflow lobe 
encasing the [Fe 11 ] jet to the east and west of the HH 666 
driving source. Proper motions reveal the dynamical age of 
these Ha outflow lobes to be ~ 1000 — 1500 yrs (see Ta¬ 
ble 2). This points to a sudden increase in the jet mass-loss 
rate ~ 1000 yrs ago that provided a kick to the ambient 
material now seen as a Hubble flow in Ha emission. Several 
more distant knots in HH 666 in dicate that the jet had o ther 
outbursts like this in the past ([Smith et alJl2003 . l2010h . 

Steady [Fe 11 ] emission inside the Ha cocoons on either 
side of the jet indicates that the high-mass-loss rate episode 
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Figure 5. Top: Colour image of HH 666 that illustrates the spatial offset between different emission lines in the jet. Bottom: Cartoon 
of the externally irradiated jet and outflow. Red lines show the collimated [Fe il] jet body that clears cavities on either side of the jet. 
Ha emission, shown in black, traces the ionized edges of these cavities and the ionized skin of the bare jet. 


has persisted, with a duration comparable to its age. A 
burst duration of ~ 1000 yrs is an order of magnitude larger 


than the estimated de 

Hartmann & Kenvon 

der of the time betw< 
by low-mass stars (e.£ 

;cay time of FU Orionis outbursts (e.g. 
1996lj. This long duration is on the or- 
een ejections found in HH jets driven 

i. Reiourth et al. 19891: Hartigan et al.l 

Il990l: Reiourth et alJ 

199_2|), and similar to the variability 


time-scale estimated for HH 46/47 ([Raga &; Noriega-Cresool 
Il992h . While the burst duration is long compared to those in¬ 
ferred for low-mass protostars, the ~ 1000 yr dynamical time 
of the inner jet is only a small fracti on of the overall age o f 
the HH 666 outflow (~ 40, 000 yr, see lReiter &; Smithll2014l h 
This points to a strong, sustained increase in the accretion 


rate from a relatively evolved intermediate-mass protostar. 
A larger sample of intermediate-mass protostars is required 
to determine the frequency and duty cycle of such accre¬ 
tion bursts, and whether they remain stronger than their 
low-mass counterparts throughout their evolution. 

[Fe n] emission from HH 666 traces a steady jet with 
higher densities and v elocities than traced by Ha emission. 
iReiter fe Smithl (1201 31 found that the mass-loss rate implied 
by the survival of Fe + in the jet is an order of magnitude 
higher than the mass-loss rate d erived from the H a emis¬ 
sion measure ([Smith et al.l l2010h . New [Fe n] spectra pre¬ 
sented here reveal velocities ±200 km s -1 throughout the 
inner jet, similar to the fastest velocities in the Ha spec- 
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Figure 6. Tracing through the HST-WFC3 Ha and [S il] images 
of HH 666. Top panel shows the Ha image of the area used to plot 
the tracings through the jet environment shown below. The plot¬ 
ted intensity is the total of the emission through the y-direction of 
the image above (perpendicular to the jet propagation direction). 
A dark cavity along the northwest limb of the jet is apparent in 
optical and IR images of the jet. This corresponds to the dip in 
the flux to the right of the driving source. The dip to the left in¬ 
dicates the presence of a complementary evacuated cavity along 
the southeast limb of the jet (see Section [472]) . 


trum. Applying the median tilt angle of the jet derived from 
the Ha velocities, a ~ 42° (see Table 2), implies that the 
three-dimensional velocity of the [Fe n] jet is ~ 300 km s _1 , 
almost twice the speed of the Ha jet. The higher mass-loss 
rate and velocities in the [Fe n] jet suggest that this compo¬ 
nent of the outflow has at least an order of magnitude more 
momentum. However, the steady velocity of the [Fe n] jet 
traveling in the wake of the initial bow shock suggests little 
momentum transfer to the local environment. This means 
that the high momentum of the [Fe n] jet will not be de¬ 
posited locally in the molecular globule, but rather injected 
into the H n region where its contribution to driving turbu¬ 
lence in the region will be dwarfed by feedback from massive 
stars. 


6 CONCLUSIONS 

In this paper, we propose a solution to the discrepancy ob¬ 
served in the morphology and kinematics of Ha and [Fe n] 
emission from HH 666, an externally irradiated protostellar 
outflow in the Carina nebula. Previous high-resolution imag¬ 
ing revealed different morphologies traced by the two lines. 
We present new, near-IR spectra of the inner jet demon¬ 
strating that [Fe n] emission also traces different kinemat¬ 
ics. Doppler velocities in the new [Fe n] spectra of the in¬ 
ner jet, features HH 666 M and O, remain remarkably con¬ 
sistent on either side of the jet, tracing velocities of ±200 


km s _1 throughout. This stands in stark contrast to the 
velocity structure seen from the sa me portio ns of the jet 
in the longslit Ha spectrum from ISmith et al.l (|2004h . Both 
HH 666 M and O show Hubble-like velocities in the Ha spec¬ 
trum with only the fastest Doppler velocities in the Hubble 
wedges on either side of the driving source approaching the 
velocities traced by [Fe n]. 

Ha and [Fe n] emission also trace different morpholo¬ 
gies, complem ent ing t he velocity discrepancies in spectra 
([Reiter fe Smithl [ 2 OI 3 I ). Sheath-like Ha emission encom¬ 
passes narrow [Fe 11 ] emission on either side of the jet. We 
argue that broad Ha emission traces the material swept up 
via prompt entrainment by a jet outburst that began ~ 1000 
yrs ago, with bright Ha arcs tracing the head of the bow 
shock that dominates entrainment. Constant velocity [Fe 11 ] 
emission, inside the outflow sheath and Hubble flow traced 
by Ha, indicates that the high-mass-loss rate episode from 
the jet is ongoing, tracing a strong accretion burst that has 
been sustained longer than the typical decay time of FU Ori- 
onis outbursts. A new epoch of Ha imaging from HST /ACS 
provides a ~ 9 yr time baseline to measure the motion of 
bright Ha features. Proper motions trace the lateral expan¬ 
sion of the Ha sheath, consistent with an expanding shell 
around the jet. 

We propose that the same underlying physics that 
drives the molecular outflows seen from more embedded 
intermediate-mass protostars creates the two-component 
outflow we observe in HH 666. Unlike jets with molecular 
outflows, HH 666 resides in a giant H II region, so both the 
jet and the material it entrains are subject to the harsh UV 
background created by the many O-type stars in the Carina 
nebula. This leads to Ha emission from the walls of the out¬ 
flow cavity and [Fe 11 ] emission from the obscured jet. Out¬ 
side the parent pillar where a cocoon no longer envelopes 
the jet, Ha and [Fe 11 ] emission trace the same morphol¬ 
ogy. Thus, HH 666 provides a missing link that unifies the 
irradiated outflows driven by unobscured protostars with 
the molecular outflows typically observed from embedded 
intermediate-mass protostars. 
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